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Abstract. We derive photometric, structural and dynamical evolution-related parameters of 11 nearby open clus¬ 
ters with ages in the range 70 Myr to 7 Gyr and masses in the range « 400 Mq to « 5 300 M©. The clusters 
are homogeneously analysed by means of J, H and Ks 2MASS photometry, which provides spatial coverage wide 
enough to properly take into account the contamination of the cluster field by Galaxy stars. Structural parameters 
such as core and limiting radii are derived from the background-subtracted radial density profiles. Luminosity 
and mass functions (MFs) are built for stars later than the turnoff and brighter than the 2MASS PSC 99.9% 
completeness limit. The total mass locked up in stars in the core and the whole cluster, as well as the correspond¬ 
ing mass densities, are calculated by taking into account the observed stars (evolved and main sequence) and 
extrapolating the MFs down to the H-burning mass limit, 0.08 M©. We illustrate the methods by analysing for 
the first time in the near-infrared the populous open clusters NGC2477 and NGG2516. For NGG2477 we derive 
an age of 1.1 ± 0.1 Gyr, distance from the Sun d© = 1.2 ±0.1 kpc, core radius Rcore = 1.4 ±0.1 pc, limiting radius 
Riim = 11.6 ± 0.7pc and total mass mtot ~ (5.3 ± 1.6) x 10® M©. Large-scale mass segregation in NGC2477 is 
reflected in the significant variation of the MF slopes in different spatial regions of the cluster, and in the large 
number-density of giant stars in the core with respect to the cluster as a whole. For NGC 2516 we derive an age 
of 160 ± 10 Myr, d© = 0.44 ± 0.02 kpc, Rcore = 0.6 ± 0.1 pc, Rum = 6.2 ± 0.2 pc and mtot (1.3 ± 0.2) x 10® M©. 
Mass-segregation in NGG2516 shows up in the MFs. Six of the 11 clusters present a slope break in the MF 
occurring at essentially the same mass as that found for the field stars in Kroupa’s universal IMF. The MF break 
is not associated to cluster mass, at least in the clusters in this paper. In two clusters the low-mass end of the 
MF occurs above the MF break. However, in three clusters the MF break does not occur, at least for the mass 
range m > 0.7 Mq. One possibility is dynamical evolution affecting the MF slope distribution. We also search for 
relations of structural and evolutionary parameters with age and Galactocentric distance. The main results for 
the present sample are: fij cluster size correlates both with age and Galactocentric distance; (ii) because of size 
and mass scaling, core and limiting radii, and core and overall mass correlate; (Hi) massive (m > lOOOM©) and 
less-massive clusters follow separate correlation paths on the plane core radius and overall mass; (iv) MF slopes 
of massive clusters are restricted to a narrow range, while those of the less-massive ones distribute over a wider 
range. Core and overall MF flattening is related to the ratio (r) of age to relaxation time. For large values of r the 
effects of large-scale mass segregation and low-mass stars evaporation can be observed in the MFs. In this sense, 
r appears to characterize the evolutionary state of the clusters. We conclude that appreciable slope flattenings in 
the overall MFs of the less-massive clusters take ~ 6 times longer to occur than in the core, while in the massive 
clusters they take a time ~ 13 times longer. We investigate cluster parameters equivalent to those determining 
the fundamental plane of ellipticals. These parameters are: overall mass, projected mass density and core radius. 
We conclude that in the present sample there is evidence of a fundamental plane. Larger samples are necessary 
to pin down this issue. 
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1. Introduction 

Open clusters are self-gravitating stellar systems which 
span a broad range of ages and masses, and present sig¬ 
nificant variations in linear dimensions and morphology 
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(e.g. Lynga II 9821 Binney & Merrifield II 9981 Tadross et 
ai.Enm. The structure of most open clusters can be de¬ 
scribed by two subsystems, the dense core and the sparse 
halo (corona). The stellar content of an open cluster pre¬ 
sumably formed at the same time from the same parent 
interstellar cloud, thus sharing similar initial conditions. 
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Stellar and dynamical evolution, as well as external inter¬ 
actions with Galactic structures, change the morphology 
and the internal mass distribution of clusters (Bergond, 
Leon & Guibert ETOTll . As a consequence of large-scale 
mass segregation high-mass stars (as well as multiple sys¬ 
tems) tend to be more concentrated in the core of evolved 
clusters while low-mass stars are transferred to the halo. 
In the oldest clusters (or the less-massive ones) the low- 
mass stars may evaporate from the cluster and get dis¬ 
persed into the surroundings (de la Fuente Marcos flHH8|l . 
Consequently, changes which occur along time must be re¬ 
flected on observational parameters such as core and over¬ 
all radii, the mass locked up in stars in the core and in the 
whole cluster, and the mass densities in both structures. 
In addition, core and overall mass-function slopes change 
as the clusters evolve (e.g. Kroupa onnri) . In this sense, 
open clusters are excellent laboratories in which models of 
star formation, stellar evolution and dynamical evolution 
can be tested both on theoretical grounds and numerical 
simulations. 

The relatively small number of stars (as compared to 
globular clusters) in open clusters poses some technical 
challenges to N-body simulations, particularly associated 
to the (i) increased granularity of the gravitational po¬ 
tential, (ii) significant fraction of binaries formed along 
the evolution, (in) significant mass loss during the course 
of the normal stellar evolution, which affects cluster life¬ 
time, and (iv) non-spherical spatial shape. A discussion 
on the above points is in e.g. de la Fuente Marcos & de 
la Fuente Marcos cnna. Gonsequently, it is important 
to search for correlations among observational parameters 
in order to better constrain theoretical studies and sim¬ 
ulation codes. Accordingly, large spatial coverage, photo¬ 
metric uniformity and a homogeneous analysis method are 
essential to obtain a consistent set of intrinsic parameters. 
This point is important since most open clusters are pro¬ 
jected against the Galactic plane and their fields may be 
heavily contaminated by Galactic (field) stars. 

Previous searches for correlations involving open clus¬ 
ter parameters can be found in, e.g. Lynga dUHUi, which 
is a compilation of literature data, Tadross et al. (Unoi, 
which is based on UBV GGD photometry, and Nilakshi, 
Pandey & Mohan II2()()2II . which is based on star counts ex¬ 
tracted from Digital Sky Survey plates. Because of the lim¬ 
ited spatial coverage and/or limited photometrical depth 
of the previous observational techniques, most parameters 
were restricted to the central regions of the clusters and 
consequently, they suffered from poor background contam¬ 
ination subtraction, as well as cluster field crowding. 

Parameters used in the present work are: (i) 
Galactocentric distance (doc)? which is an external pa¬ 
rameter; (ii) cluster age, which is associated to the stel¬ 
lar evolution; (Hi) core (Rcore) and limiting (Rum) radii, 
which are structural parameters associated to the dynam¬ 
ical evolution; (iv) core and overall mass; (v) core and 
overall projected mass density (in Mq pc“^); (vi) core and 
overall mass density (in Mqpc”^). Parameters listed in 
(ii) and (iv) to (vi) are intrinsically associated to the stel¬ 


lar and dynamical evolution. We refer to the whole struc¬ 
ture (0 < r < Riini) of a cluster as overall. 

For the sake of spatial and photometric uniformity, we 
employ in the present work J, H and Ks 2MASS^ photom¬ 
etry. The 2MASS Point Source Catalogue (PSC) is uni¬ 
form reaching relatively faint magnitudes covering nearly 
all the sky, allowing a proper background definition even 
for clusters with large angular sizes. The clusters in the 
present work have already been studied by us in previous 
papers (in proper motion and dynamical state studies), ex¬ 
cept for the rich open clusters NGC2477 and NGC2516, 
which are studied in the near-infrared for the first time and 
illustrate in detail the present analysis method as well. 

Age and Galactocentric distance (derived from the dis¬ 
tance from the Sun) are obtained by fitting isochrones 
to the near-infrared 2MASS colour-magnitude diagrams 
(CMDs), which gives as well the foreground colour ex¬ 
cess in that direction (e.g. Bica, Bonatto & Dutra 120133 
I2flfl4hll . The statistical significance of these parameters de¬ 
pends directly on the quality and depth of the photometry 
(Bonatto, Bica & Pavani E()()4l Bonatto, Bica & Santos Jr. 
i™il . The core radius is derived by fitting a King model 
to the radial density profile, while the limiting radius cor¬ 
responds to the distance from the center where the radial 
density profile merges with the background level. In this 
case, deep enough photometry is necessary to produce a 
radial density profile which includes low main-sequence 
(MS) stars as well, while spatial coverage is fundamen¬ 
tal, since a high background may produce an underesti¬ 
mated cluster dimension (Bica, Bonatto & Dutra l^04bl 
Bonatto, Bica & Pavani 1^0041 Bonatto, Bica & Santos Jr. 
120011 ). Wide spatial coverage is important for the dynam¬ 
ically evolved clusters, since a significant fraction of the 
low-mass content may already have been dispersed into 
the surroundings (Bonatto & Bica 120131 Bonatto, Bica & 
Pavani 1^0041 Bonatto, Bica & Santos Jr. 1200411 . Core and 
overall masses are obtained by first deriving the intrin¬ 
sic luminosity (LF) and mass (MF) functions, which de¬ 
pend essentially on properly taking into account the stel¬ 
lar background distribution, and then integrating the MFs 
down to the H-burning limit stars. Consequently, depth 
and coverage are fundamental to properly derive mass, 
particularly for clusters projected against the disk and/or 
bulge (Bica, Bonatto & Dutra I2004bl Bonatto, Bica & 
Santos Jr. l2004|l . 

In the present work we intend to search for funda¬ 
mental relations in the parameter space described above. 
Additionally, since star clusters are, to a first approxi¬ 
mation, scaled-down versions of elliptical galaxies, we in¬ 
tend to explore the possibility of an analogous fundamen¬ 
tal plane (FP, e.g. Dressier et al. llH87l Djorgovski & Davis 
ITWIi of open cluster parameters. The parameters which 
could correspond to those of the FP in ellipticals could be 
cluster overall and/or core mass, core radius (the analo- 

^ The Two Micron All Sky Survey, All Sky 
data release (Skrutskie et al. I1997II . available at 
http://www.ipac.caltech.edu/2mass/releases/allsky/ 
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gous of the effective radius in the ellipticals), and the pro¬ 
jected mass density. However, the conceptual differences 
between ellipticals and open clusters, e.g. total popula¬ 
tion and dynamical evolutionary state, may hinder such 
eventual relations. 

This paper is organized as follows. In Sect. |2 we pro¬ 
vide general data on the target clusters. In Sect. 0 we illus¬ 
trate the analysis methods with NGC 2477 and NGC 2516. 
Gluster parameters are presented and correlations are dis¬ 
cussed in Sect. 01 The MF slope break and dynamical ef¬ 
fects on MFs are discussed in Sect.|Sl In Sect.|Slwe discuss 
a possible FP in open clusters. Finally, concluding remarks 
are given in Sect. 0 

2. The target clusters 

We selected for the present work 9 open clusters already 
studied by us with 2MASS in previous papers for differ¬ 
ent goals: M 26 (NGG 6694), NGC 2287, M 48 (NGC 2548), 
M93 (NGC 2447), NGC 5822, NGC 3680, IC4651, M67 
(NGC 2682) and NGC 188. Proper motions and binarity 
in M26, NGC 2287, M48, M93, NGC 5822, NGC 3680, 
IC4651 and M67 were investigated in Bica & Bonatto 
cnni. The advanced dynamical state of NGC 3680 was 
analysed in Bonatto, Bica & Pavani (mni . Mass segrega¬ 
tion in M67 was studied in Bonatto & Bica (j2IS|). The 
dynamical state of the old open cluster NGC 188 was in¬ 
vestigated in Bonatto, Bica & Santos Jr. 

The populous open clusters NGC 2477 and NGC 2516 
are discussed in Sect. 0 The clusters above span a wide 
age range, which is important since age is a fundamental 
factor associated to dynamical evolution. General data for 
the above open clusters are given in Table 0 The data 
in cols. 2-6 are from the WEBDA^ open cluster database 
('Mermilliod ll 996]l . In cols. 7-12 we list parameters derived 
from the 2MASS photometry. 

We indicate below papers and results relevant to the 
present work, discriminated by open cluster. 

M 26: Battinelli, Brandimarti & Capuzzo-Dolcetta 
give a limiting radius of Rum « 3.3 pc and a total stellar 
mass of mtot ~ 240 Mq. 

NGC 2516: Battinelli, Brandimarti & Capuzzo-Dolcetta 
give Riim ~ 1.9 pc and mtot ~ 282 Mq. Bergond, 
Leon & Guibert (mnn) give d© « 0.4 kpc, an age of 
« 113 Myr, Rtidai = 1.9pc and mtot ~ 170 Mq. Tadross et 
al. llWa give d© « 0.4 kpc, age « 63 Myr, Rum « 1.8 pc, 
dec ~ 8.5kpc and mtot ~ 36 M©. Jeffries, Thurston & 
Hambly Ij2001ll estimate an age of « 150 Myr and total 
mass of mtot ~ 1240—1560 Mq. They also detected a sharp 
change in the MF slope for stars with 0.3 < m(MQ) < 0.7. 

NGC 2287: Battinelli, Brandimarti & Capuzzo-Dolcetta 
give Riim « 4.4 pc and mtot ~ 195 Mq. Bergond, 

^ http://obswww.unige.ch/webda 


Leon & Guibert llWTll included this object in a study 
of the gravitational tidal effects on open clusters. They 
give a distance from the Sun of dQ ~ 0.7kpc, an age of 
« 240 Myr, a tidal radius of Rtidai = 4.1 pc and a total 
stellar mass of mtot ~ 300 Mq. 

M 48: Bergond, Leon & Guibert J^nTTlTi give dQ « 0.8 kpc, 
an age of « 360 Myr and Rtidai = 4.8 pc. 

M 93: Nilakshi, Pandey & Mohan II2002II derived Rcore ~ 
0.5 pc and Rum ~ 3.3 pc. 

NGC 5822: Battinelli, Brandimarti & Capuzzo-Dolcetta 
mm give mtot « 251 Mq. 

NGC 2477: Based on UBV photometry of 83 red-giant 
candidates in the field of NGC 2477, Eigenbrod et al. 
TO derived a core radius Rcore = 1.8 pc and a tidal 
radius Rtidai = 8.1 pc. They also estimate that the total 
mass locked up in stars more massive than 1 Mq must 
be larger than 5 400 Mq . A break in the mass function in 
the lower main sequence at m « 0.27 Mq was suggested 
by von Hippel et al. Nilakshi, Pandey & Mohan 

coni derived Rcore ~ 2.2 pc and Rum ~ 9.1 pc. Tadross 
et al. II2()02II givedQ « 1.4 kpc, age « 1 Gyr, Rum « 5.5 pc, 
dec ~ 9.0kpc and mtot ~ 2100 Mq. 

NGC 3680: Bonatto, Bica & Pavani mm studied this dy¬ 
namically advanced open cluster and derived dQ « 1 kpc, 
age « 1.6Gyr, Rcore « 0.7pc, Rum « 6.4pc and mtot « 
550 Mq. 

IC 4651: Meibom, Andersen & Nordstrom mm derived 
an age of ~ 1.7 Gyr, Rtidai ~ 5.8 pc and a present total 
mass of ~ 630 Mq; they also present evidence of moderate 
mass segregation. Tadross et al. II2002II give dQ « 0.8 kpc, 
age « 2.2 Gyr, Rum « 1.4 pc, dec « 7.8 kpc and mtot « 
100 Mq. 

M 67: Battinelli, Brandimarti & Capuzzo-Dolcetta mm 
give mtot ~ 407 Mq. Bonatto & Bica mm derived an 
age of « 3.2 Gyr, dQ « 0.9 kpc, Rcore = 1.2 pc and 
Riim ~ 1.8 pc. They also detected mass segregation. Mass 
segregation in M67 was also detected by Sarajedini et 
al. II1999II . Montgomery, Marschall & Janes II1993II esti¬ 
mate a mass of « 800 Mq. Detailed N-body modelling 
of the dynamical and stellar evolution of M67 was car¬ 
ried out by Hurley et al. mm : they estimate a much 
larger initial mass for M67. Tadross et al. mm give 
dQ « 0.8 kpc, age « 4 Gyr, Rum ~ 3.2 pc, dec ~ 9.1 kpc 
and mtot ~ 820 Mq. 

NGC 188 : Bonatto, Bica & Santos Jr. mmi derived an 
age of « 7Gyr, dg « 1.7kpc, Rcore ~ 1.3pc, Rtidai ~ 
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Table 1. General data on the target clusters. 


Cluster 

(1) 


WEBDA 



Present work 




a(2000) 

(hms) 

(2) 

J(2000) 

(“'") 

(3) 

Age E(B-V) do 
(Myr) (kpc) 

(4) (5) (6) 

a(2000) 5(2000) 

(hms) (°'") 

(7) (8) 

Age 

(Myr) 

(9) 

E(B - V) 

(10) 

d© 

(kpc) 

(11) 

dec 

(kpc) 

(12) 

M26 

18 : 45 : 18 

-09 : 23 : 00 

85 

0.59 

1.60 

18 : 45 : 14^ -09 : 22 : 55’' 

70± 10 

0.42 ±0.03 

1.6 ±0.1 

6.6 

NGC 2516 07 : 58 : 04 

-60 : 45 : 12 

113 

0.10 

0.41 

07 : 58 : 04 -60 : 45 : 12 

160 ± 10 

0.00 

0.44 ± 0.02 

8.0 

NGC 2287 06:46: 01 

-20 : 45 : 24 

243 

0.03 

0.69 

06 : 46 : 01 -20 : 45 : 24 

160 ± 10 

0.00 

0.8 ±0.1 

8.5 

M48 

08 : 13 : 43 

-05 : 45 : 00 

360 

0.03 

0.77 

08 : 13 : 45^ -05 : 46 : 42^ 

360 ± 40 

0.00 

0.8 ±0.1 

8.5 

M93 

07 : 44 : 30 

-23 : 51 : 24 

387 

0.05 

1.04 

07 : 44 : 31^ -23 : 50 : 42^ 

400 ± 50 

0.03 ±0.01 

1.1 ±0.1 

8.6 

NGC 5822 15 : 04 : 21 

-54 : 23 : 48 

662 

0.15 

0.92 

15 : 04 : 21 -54 : 23 : 48 

1000 ± 100 

0.00 

0.7 ±0.1 

7.4 

NGC 2477 07 : 52 : 10 

-38 : 31 : 48 

705 

0.28 

1.22 

07 : 52 : 18'^ -38 : 31 : 57^ 

1100 ± 100 

0.06 ± 0.03 

1.2 ±0.1 

8.4 

NGC 3680 11 : 25 : 38 

-43 : 14 : 36 

1200 

0.07 

0.94 

11 : 25 : 38 -43 : 14 : 36 

1600 ± 100 

0.00 

1.0 ±0.1 

7.8 

IC4651 

17 : 24 : 49 

-49 : 56 : 00 

1140 

0.12 

0.89 

17 : 24 : 521‘ -49 : 56 : 52^ 

1800 ± 200 

0.00 

0.9 ±0.1 

7.2 

M67 

08 : 51 : 18 

-kll : 48 : 00 

2560 

0.06 

0.91 

08 : 51 : 16^ -til : 48 : 54^ 3200 ± 100 

0.00 

0.9 ±0.1 

8.7 

NGC 188 

00 : 47 : 28^ -k85 : 15 : ISt 

4300 

0.08 

2.05 

00 : 47 : 53 -t85 : 15 : 30 7000 ± 1000 

0.00 

1.7 ±0.1 

8.9 


Table Notes. Colour excesses in col. 10 are derived from the isochrone fit to the 2MASS CMDs. Uncertainty in Galactocentric 
distance in col. 12 is « 0.1 kpc for all clusters. Uncertainties in cols. 9 to 12 are derived from the isochrone fit. (^) Coordinates 
optimized to maximize the number-density of stars in the cluster center (Sect.|^. 


21 pc and a total mass of « 4 000 Mq; they detected strong 
mass segregation. Mass segregation was also detected by 
Sarajedini et al. II1999II . Nilakshi, Pandey & Mohan 1200211 
derived Rcore ~ 2.1 pc and Rum « 6.5 pc. Tadross et al. 
1200211 give d© = 1.5 kpc, age « 5Gyr, Rum « 3.1 pc, 
dec ~ 9.4kpc and an exceedingly low mass of « 208 Mq. 

The near-infrared GMDs and radial density profiles 
of M26, NGG2287, M48, M93, NGC5822, NGC3680, 
1C 4651 and M 67 are shown in Bica & Bonatto (EnnH). For 
NGC 188 they are in Bonatto, Bica & Santos Jr. 1200411 . 

Comparing WEBDA compilation with the present re¬ 
sults (Table P), ages and distances from the Sun agree in 
general. However, the colour excesses E(B — V) tend to be 
smaller when derived from the near-infrared, as discussed 
in Bica & Bonatto 1200411 . 


3. The populous open clusters NGC 2477 and 
NGC 2516 

3.1. The intermediate-age open cluster NGC2477 

According to WEBDA, the central coordinates of 
NGC 2477 are (J2000) a = 07'‘52™10^ and S = 
—38°31'48". However, the corresponding radial density 
profile (built with a step in radius of Ar = 1') presented a 
dip for r = Ob We searched for a new center by examining 
histograms for the number of stars in bins of right as¬ 
cension and declination. The resulting coordinates which 
maximize the density of stars at the center are (J2000) 
a = 07^52^18", and 6 = -38°31'57", corresponding to 
£ = 253.58° and b = —5.82°. In what follows we refer to 
these coordinates as the center of NGC 2477. 


3.1.1. The 2MASS photometry and near-infrared CMD 

The VizieR^ tool has been used to extract 2MASS pho¬ 
tometry of the stars in a circular area with radius r = 60' 
centered on the optimized coordinates (Tabled]). 2MASS 
photometric uncertainties as a function of magnitude are 
discussed in Soares & Bica and Bonatto, Bica, & 

Santos Jr. 1200411 . 

Because of the moderately high Galactic latitude {b « 
—6°), the field of NGC 2477 presents limited background 
contamination. This can be seen in Fig. ^ in which we 
show the J X (J — H) CMD of the central (r = 10' ~ 3.6 pc) 
region of NGC 2477 (left panel) compared with the corre¬ 
sponding, same area, offset field (right panel). To match 
the areas of both CMDs in Fig. ^the offset field has been 
obtained by extracting stars in the outermost ring, lo¬ 
cated at 59.16' < r < 60', which represents essentially the 
Galactic background/foreground stellar contribution. The 
CMD of the central region of NGC 2477 presents a well- 
defined morphology including a densely populated MS up 
to the turnoff and a rich giant clump that strongly con¬ 
trasts with the CMD of the offset field. We show in both 
panels the colour-magnitude filter which is used to in¬ 
crease membership probability by discarding stars with 
far-off colours with respect to the isochrone fit (see be¬ 
low). 

Cluster parameters have been derived by fitting solar 
metallicity Padova isochrones (Girardi et al. 1200211 com¬ 
puted with the 2MASS J, H and Ks filters'^, to the ob¬ 
served CMD in Fig. ^ The 2MASS transmission filters 
produced isochrones very similar to the Johnson ones, 
with differences of at most 0.01 in (J — H) (Bonatto, Bica 
& Girardi I2004II . For reddening and absorption trans¬ 
formations we use Rv = 3.2, and the relations Aj = 


® http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=n/246 

4 
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(J-H) (J-H) 


Fig. 1. Left panel: J x (J — H) CMD of the central 10' 
(~ 3.6 pc) of NGC 2477 with the best-fit Padova isochrone 
(age ~ 1.1 Gyr) superimposed. Note the densely popu¬ 
lated main sequence and rich giant clump. Right panel: 
CMD of the corresponding (same area) offset field. The 
colour-magnitude filter is shown in both panels. 

0.276x Ay and E(J — H) = 0.33 x E(B — V), according to 
Dutra, Santiago & Bica and references therein. The 

best-fit for NGC 2477 uses the 1.12 Gyr isochrone. Taking 
into account the uncertainties associated to the isochrone 
fit we derive for NGC 2477 an age of 1.1 ± 0.1 Gyr, 
a colour excess E(J — H) = 0.02 ± 0.01 which corre¬ 
sponds to E(B — V) = 0.06 ± 0.03, a distance modulus 
(m — M)j = 10.5 ± 0.1, and a distance from the Sun 
d© = 1.23 ± 0.06 kpc. With this value the Galactocentric 
distance of NGC 2477 turns out to be dec = 8.4±0.1kpc, 
using 8.0 kpc as the distance of the Sun to the center of 
the Galaxy (ReidESHl- 

3.1.2. Structural parameters 

Structural parameters of NGC 2477 have been derived by 
means of the radial density profile, defined as the pro¬ 
jected number of stars per area around the cluster center. 
Before counting stars we applied the colour filter shown in 
Fig. [2 to the CMD of the cluster (for stars with distance 
to the center from 0' to 60') in order to discard most of 


the contamination by Milky Way stars. This procedure 
has been previously applied in the analysis of the open 
clusters M67 (Bonatto & Bica 1200,^ . NGC 188 (Bonatto, 
Bica & Santos Jr. I2004|l and NGC 3680 (Bonatto, Bica 
& Pavani EOOll). To further minimize the probability of 
background contamination, spurious detections and the 
increase of photometric uncertainties at faint magnitudes 
we restricted the radial density analysis to filtered stars 
with J < 14.5, corresponding to Mj < 4.0 and mass 
m > 0.95Mq ISect. ITl .3|l . At J = 14.5 the photomet¬ 
ric uncertainties amount to ej ~ 0.04 (Bonatto, Bica & 
Santos Jr. Enoni. The radial density profile was obtained 
by counting stars inside concentric annuli with a step of 
1.0' in radius. The background contribution level corre¬ 
sponds to the average number of stars included in the 
ring located at 40' < r < 60' (14.3 < r(pc) < 21.5). The 
above procedures were applied to the filtered CMD. 

The resulting radial density profile of NGC 2477 is 
shown in Fig. [21 For absolute comparison between clus¬ 
ters we scale the radius in the abscissa in parsecs, and the 
number density of stars in the ordinate in stars pc“^ using 
the distance derived in Sect. uni The statistical signifi¬ 
cance of the profile is reflected in the relatively small Icr 
Poisson error bars. 

Structural parameters of NGC 2477 were derived by 
fitting the two-parameter King Ill966a[l surface density 
profile to the background-subtracted radial distribution of 
stars. The two-parameter King model essentially describes 
the intermediate and central regions of normal clusters 
CKing fr^bbl Trager, King & Diorgovski ll995|l . The fit was 
performed using a nonlinear least-squares fit routine which 
uses the errors as weights. The best-fit solution is shown 
in Fig. 121 superimposed on the radial density profile. From 
the fit we derived a core radius Rcore = 1.4 ± 0.1 pc and a 
projected central density of stars of cto « 70± 7 stars pc“^. 
Our value of Rcore is somewhat smaller than Rcore = 1-8 pc 
of Eigenbrod et al. (Enni) and Rcore = 2.2 pc of Nilakshi, 
Pandey & Mohan II2()()2II . The present central density of 
stars is ~ 3x that of Nilakshi, Pandey & Mohan mm . 

Within uncertainties the two-parameter King model 
reproduces well the radial density profile in NGC 2477. 
Since it follows from an isothermal (virialized) sphere, 
the close similarity of the radial distribution of stars in 
NGC 2477 with a King profile suggests that despite its in¬ 
termediate age, this cluster has evolved into an advanced 
dynamical state. We will return to this point in Sect. ET21 

Considering the radial density profile fluctuations with 
respect to the background level, we can define a limiting 
radius (Rum) for the cluster, in the sense that for regions 
beyond Riim, the null-contrast between cluster and back¬ 
ground star density would produce large Poisson errors 
and consequently, meaningless results. Thus, for practi¬ 
cal purposes, the bulk of the cluster’s stars are contained 
within Riini. In the case of NGC 2477 the limiting radius 
turns out to be Rum = 11.6±0.7pc. Finally, we derived the 
tidal radius (Rtidai) by fitting the three-parameter King 
(CMl model to the background-subtracted radial distri¬ 
bution of stars in Fig. (21 The resulting Rtidai = 16 ± 2 pc 
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Fig. 2. Projected radial density profile of stars in 
NGC 2477. The average background level is shown as the 
narrow shaded rectangle; Icr Poisson errors are also shown. 
The solid line shows the two-parameter King model fit to 
the radial profile. The extent of the core region and that 
of the offset field are indicated. 

is ~ 33% larger than Riim- Our value of Rtidai is ~ 2 times 
larger than the values given by Eigenbrod et al. mm 
and Nilakshi, Pandey & Mohan 1121)0211 . and a factor of 
^ 3 larger than that of Tadross et al. 1I2()()2II . The sig¬ 
nificant differences between our values of Rcore and Rtidai 
and those of previous works may be accounted for by the 
depth of the 2MASS photometry, the proper background 
definition and the use of a colour-magnitude filter. 

3.1.3. Mass functions 

The populous nature of NGC 2477 provides an opportu¬ 
nity to study in detail with a high statistical significance 
the spatial distribution of mass functions = ^). 

Based on the King profile fit and taking into account the 
features present in the radial density proHle of NGC 2477 
(Fig.ia we selected the following regions: (i) 0.0 < r(pc) < 
1.4 (core), (ii) 1.4 < r(pc) < 5.2, (Hi) 5.2 < r(pc) < 11.6 
and (iv) 0.0 < r(pc) < 11.6 (overall). Regions (ii) and (Hi) 
are the halo of NGC 2477. In order to maximize the sig¬ 
nificance of background counts, we consider as offset held 
the outermost ring at 14.3 pc < r < 21.5 pc, which lies 
^ 2.7pc beyond the cluster’s limiting radius. 


Although relatively small, the Galactic stars contam¬ 
ination of the CMD of NGC 2477 must be taken into ac¬ 
count in order to derive the intrinsic luminosity and mass 
distributions of the member stars. To do this we hrst ap¬ 
ply the colour hlter (Fig.^J to both cluster and offset held 
CMDs. The hltering process takes into account most of 
the background, leaving a residual contamination. We deal 
with this residual contamination statistically by building 
the LFs for each cluster region and offset held. We build 
LFs for the three 2MASS bands independently, taking into 
account the 99.9% Point Source Catalogue Completeness 
Limit^. Consequently, the faint magnitude limit of each 
LF is J = 15.8, H = 15.1 and Ks = 14.3, respectively. We 
take the turnoh as the bright limit to avoid inconsistencies 
in the mass-luminosity relation. For each 2MASS band we 
build a LF by counting stars in magnitude bins from the 
respective faint magnitude limit to the turnoff, both for 
each cluster region and offset held. Considering that the 
solid angle of the offset held is diherent from that of a given 
cluster region, we multiply the offset held LF by a numer¬ 
ical factor so that the solid angles match. The intrinsic 
LF of each cluster region is obtained by subtracting the 
respective (i.e. solid angle-corrected) offset held LF from 
that of the cluster region. Finally, the intrinsic LFs are 
transformed into MFs using the mass-luminosity relation 
obtained from the 1.12 Gyr Padova isochrone and the dis¬ 
tance modulus (m — M)j = 10.5. Remark that these pro¬ 
cedures are repeated independently for the three 2MASS 
bands. The hnal MF of a given cluster region is produced 
by combining the J, H and Ks MFs into a single MF. 
The resulting spatial MFs of NGC 2477, covering the mass 
range 0.76 < m(MQ) < 1.95, are shown in Fig.|3 

The MFs present a break followed by slope battening 
for masses in the range 0.76 < m(M0) < 1.08, which is 
noticeable particularly in the core (panel (d)) and in the 
region 1.4 < r(pc) < 5.2 (panel (c)). We note that von 
Hippel et al. (dS, based on deep V and I HST and 
WFPC2 photometry of the central parts of NGC 2477, 
described a break in the MF, but occurring at a lower 
mass value, m « 0.27 Mq (Sect. El). 

Bearing in mind the MF break at mbreak = 1.08 Mq 
we ht the function (j)(m) oc in the mass ranges 

0.76 < m(Mo) < 1.08 and 1.08 < m(M0) < 1.95. The 
resulting hts are shown in Fig. El and the MF slopes of 
each mass range are given in cols. 4 and 5 of Table O 
respectively. 

We provide in Table El parameters derived from the 
LFs and MFs. The number of evolved stars (col. 2) is cal¬ 
culated by integrating the intrinsic LFs for magnitudes 
brighter than the turnoff. Multiplying this by the mass at 
the turnoff (1.98 M0) gives an estimate of the evolved- 
star mass (col. 3). This procedure produces a realistic 
value of the number of member evolved stars because the 
background contamination has already been statistically 
subtracted from the LF. The MF slopes in the two mass 

® Corresponding to the Level 1 Requirement, according to 

http ://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6^5al.html 
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Fig. 3. Mass functions in different spatial regions of 
NGC2477. Each panel contains MFs derived from the J, 
H and Ks 2MASS photometry. A break in the MFs at 
m « 1.08 Mq can be seen particularly in panels (c) and 
(d). ME fits ((/'(to) oc to“*^^+^)) for each mass range are 
shown as dashed lines. The respective MF slopes are given. 
Note that the flattest (panel (d)) and steepest (panel (b)) 
MFs coincide spatially in both mass ranges. 

ranges are in cols. 4 and 5. The number of MS stars and 
corresponding stellar mass are derived by integrating the 
MF from the faint magnitude limit to the turnoff. We add 
to these the corresponding values of the number and mass 
of evolved stars to derive the total number of observed 
stars (col. 6), observed mass (col. 7), projected mass den¬ 
sity (col. 8) and mass density (col. 9). 

An estimate of the total mass locked up in stars in 
NGC 2477 was made by taking into account all stars from 
the turnoff down to the H-burning mass limit, 0.08 Mq. 
We do this by directly extrapolating the low-mass MFs 
(0.76 < m(M0) < 1.08) down to 0.08 Mq, except for 
the region 5.2 < r(pc) < 11.6 for which the slope is 
X = 1.9±0.7. In this case we assume the universal 
Initial Mass Function (IMF) of Kroupa l|2001|l . in which 
X = 0.3 ± 0.5 for the range 0.08 < m(M0) < 0.5. In the 
range 0.5 < m(M0) < 0.76 we use our value of x which, 
within the uncertainty is similar to the value presented 
in Kroupa II2()()1II for the equivalent mass range. The re¬ 
sulting extrapolated values of the number of stars and ex¬ 


trapolated mass (added to the corresponding values of the 
number and mass of evolved stars) are given respectively 
in cols. 10 and 11 of Table [3 The extrapolated projected 
and mass densities are in cols. 12 and 13. 

The total mass locked up in stars in the core of 
NGC 2477 results in mtot(core) « (5.2 ± 0.1) x 10^ Mq, 
nearly 25% larger than the mass corresponding to the 
observed stars. The total overall mass is mtot (overall) « 
(5.3 ±1.5) X 10^ Mq, which is in close agreement with the 
mass estimate (for stars with m > IMq) of Eigenbrod 
et al. II2()()4II and ~ 2.5 times larger than the estimate 
of Tadross et al. II2()()2II . The extrapolated overall mass is 
nearly twice the mass stored in the observed stars. 

3.1.4. Mass segregation 

Within the uncertainty, the overall MF slope (y = 1.5 ± 
0.1) in the mass range 1.08 < m(MQ) < 1.95 is similar 
to a standard Salpeter (y = 1.35) IMF. However, the MF 
slope presents large variations in the inner regions, being 
flat (x = 0.0 ± 0.2) in the core and steep (x = 3.3 ± 0.24) 
in the region 5.2 < r(pc) < 11.6. MF slope flattening from 
the outskirts to the core is detected in the mass range 
0.76 < m(MQ) < 1.08 as well. Similarly to the mass range 
1.08 < m(M0) < 1.95 the flattest MF (y = -2.3 ± 0.5) 
occurs in the core, while the steepest MF (y = 1.9 ± 0.7) 
occurs in the region 5.2 < r(pc) < 11.6. This fact reflects 
the advanced dynamical state of NGC 2477, particularly 
the effects of large-scale mass segregation, in the sense that 
low-mass stars originally in the core are being transferred 
to the cluster’s outskirts while massive stars accumulate 
in the core. This produces a flat MF in the core and a 
steep one in the halo. 

Another evidence pointing to large-scale mass segre¬ 
gation is the large difference in the number-density of gi¬ 
ants in each cluster region with respect to that of MS 
stars. We assume here as giants the stars brighter than 
the turnoff (with mass m > 1.98 Mq) and as MS stars 
those with mass in the range 0.08 < m(M0) < 1.98. For 
the giants the number-density ratio is p«(dv°raii) 
while for the MS stars the ratio drops to G ^ 30. 

^ p*(overalij 

This means that in the core of NGC 2477 the number- 
density of giants with respect to MS stars is ~ 4 larger 
than in the cluster as a whole. The situation gets reversed 
in the region 5.2 < r(pc) < 11.6, where the giants and MS 
number-density ratios are, respectively ^ 0.15 

and ^ 0.61. Consequently, in this region the number- 
density of MS stars with respect to the giants is ^ 4 larger 
than in the cluster as a whole. 

Mass segregation in a star cluster scales with the 
relaxation time, defined as treiax = g^^tcross, where 
tcross = R./(rv is the crossing time, N is the (total) num¬ 
ber of stars and Uv is the velocity dispersion (Binney 
& Tremaine EMI). treiax is the characteristic time scale 
in which the cluster reaches some level of kinetic energy 
equipartition with massive stars sinking to the core and 
low-mass stars being transferred to the halo. Assuming a 
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Table 2. Data derived from the mass functions in NGC2477 and NGC2516. 








NGC 2477 







Evolved 

X 



Observed -|- Evolved 


Extrapolated -|- Evolved 


Region 

N^ 

ITlevol 

0.76 - 1.08 

1.08 - 1.95 

N* 

mobs ^ 

P 

N* 

mtot 

(7 

P 

(pc) 

(Stars) 

(lO^'Mo) 



(lO^stars) 

(lO^Mo) 


(lO^stars) 

(io^Mq) 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) (8) 

(9) 

(10) 

(11) 

(12) 

(13) 

0.0 - 1.4 

30 ± 6 

0.6 ± 0.1 

-2.3 ± 0.5 

0.0 ± 0.2 

0.3 ± 0.03 

0.4 ± 0.04 66 ±7 

37 ± 4 

0.7 ±0.3 

0.5 ± 0.1 

85 ± 13 

45 ± 7 

1.4 - 5.2 

81 ± 8 

1.6 ± 0.2 

-0.6 ± 0.4 

1.4 ± 0.2 

1.2 ± 0.07 

1.5 ±0.1 20 ±1 

2.8 ± 0.2 

5.4 ±3.3 

2.7 ± 0.6 

34 ± 8 

4.6 ± 1.1 

5.2 - 11.6 

19 ± 9 

0.4 ± 0.2 

1.9 ± 0.7 

3.3 ± 0.4 

0.7 ± 0.07 

0.7 ±0.1 2.1 ±0.2 

0.13 ± 0.01 

6.2 ±4.4 

2.2 ± 0.9 

7±3 

0.4 ± 0.1 

0.0 - 11.6 

128 ± 12 

2.5 ± 0.2 

0.4 ±0.4 

1.5 ± 0.1 

2.2 ± 0.11 

2.7 ±0.1 6.3 ±0.3 

0.43 ± 0.02 

12.1 ± 7.7 

5.3 ± 1.5 

13 ± 4 

0.8 ± 0.2 







NGC 2516 







Evolved 

X 



Observed -|- Evolved 


Extrapolated ± Evolved 


Region 

N* 

^evol 

0.39 - 0.90 

0.90 - 4.17 

N* 

niobs ^ 

P 

N* 

mtot 

(7 

P 

(pc) 

(Stars) 

(lO^'Mo) 



(lO^stars) 

(lO^'Mo) 


(lO^stars) 

(lO^'Mo) 



0.0 - 0.6 

— 

— 

-1.1 ± 0.5 

0.7 ±0.2 

0.6 ± 0.1 

0.9 ±0.1 79 ±14 

93 ± 16 

0.8 ± 0.2 

0.9 ± 0.2 

82 ± 14 

103 ± 17 

0.6 - 1.7 

— 

— 

0.3 ± 0.3 

0.7± 0.2 

1.6 ± 0.2 

1.9 ±0.3 26 ±4 

11 ± 2 

3.9 ± 1.2 

2.4 ± 0.4 

32 ± 5 

13.3 ± 2.1 

1.7 - 6.2 

— 

— 

0.5 ± 0.3 

1.8 ± 0.1 

7.4 ± 1.2 

6.4 ±0.8 5.5 ±0.6 

0.6 ± 0.1 

26.9 ± 10.3 

9.8 ± 1.8 

8.6 ± 1.6 

1.0 ± 0.2 

0.0 - 6.2 

— 

— 

0.4 ±0.2 

1.4 ±0.1 

9.6 ± 1.1 

9.1 ±0.8 7.4 ±0.7 

0.9 ± 0.1 

30.1 ± 8.7 

12.8 ± 1.6 

10.5 ± 1.3 

1.3 ± 0.2 


Table Notes. Cols. 4 and 5 give the MF slopes derived for the low-mass and high-mass MS ranges. Notice that the low-mass MS ranges are 
0.76 < m(MQ) < 1.08 and 0.39 < m(MQ) < 90, respectively for NGC2477 and NGC2416. The high-mass MS ranges are 1.08 < m(MQ) < 1.50 
and 0.90 < m(MQ) < 4.17. The mass of the evolved stars is included in mobs (col. 7) and mtot (col. 11). Units of a (cols. 8 and 12) and p 
(cols. 9 and 13) are Mq pc“^ and Mq pc~^, respectively. 


typical (Tv ~ 3kms“^ (Binney & Merrifield 1 199811 we ob¬ 
tain for the whole cluster treiax ~ 600 Myr, and for the 
core treiax ~ 6 Myr. Consequently, the presence of mass 
segregation and thus MF slope flattening, particularly in 
the core, is consistent with both estimates of treiax which 
are smaller than the age of NGC 2477. 

3.2. The moderately young open cluster NGC2516 

The analysis of NGC 2516 followed the same steps as that 
of NGC 2477, noting that in the case of NGC 2516 the 
WEBDA coordinates (Table ^1 corresponded to the ac¬ 
tual cluster center. The Galactic coordinates of NGC 2516 
are £ = 273.82° and b = —15.86°. Because of the large 
projected area of NGC 2516 we used a radius of r = 80' 
to extract the 2MASS photometry. The J x (J — H) CMD 
of the central (r = 20' ~ 2.5 pc) region of NGC 2516 is 
shown in the left panel of Fig. 0] and the corresponding, 
same area (extracted from 77.46' < r < 80') offset field is 
in the right panel. The young age of NGC 2516 is appar¬ 
ent in the CMD morphology, particularly the extended, 
nearly vertical MS and absence of evolved stars. 

3.2.1. Structural parameters 

The best-fit to the CMD of NGC 2516 was obtained with 
the 160 Myr Padova isochrone, a colour excess E(J — H) = 
0.0 (E(B — V) = 0.0), a distance modulus (m — M)j = 
8.2±0.1, and a distance from the Sun d© = 0.44±0.02kpc. 
With this value the Galactocentric distance of NGC 2516 
turns out to be dec = 8.0 ± 0.1 kpc. Considering the 
uncertainties associated to the isochrone fit the age of 
NGC 2516 gets constrained to the range 160±10Myr. This 


value agrees with that estimated by Jeffries, Thurston & 
Hambly mm . 

The radial density profile of stars in NGC 2516 (after 
applying the colour-magnitude filter - Fig. El) is shown in 
Fig. El where the background contribution corresponds to 
the average density of stars in the region 7.6 < r(pc) < 
10.18. Similarly to NGC 2477 we restricted the radial den¬ 
sity analysis to Altered stars with J < 13.0, correspond¬ 
ing to Mj < 4.8 and mass m > 0.8 Mq (Sect. 13.2.21) . 
At J = 13.0 the photometric uncertainties amount to 
ej ~ 0.03 (Bonatto, Bica & Santos Jr. 120041) . The fit with 
the two-parameter King Ill966al) surface density profile to 
the background-subtracted radial density profile resulted 
in Rcore = 0.6 ± 0.1 pc and a projected central density of 
stars of (To « 67±9 stars pc“^. We derived a limiting radius 
of Riim = 6.2 ± 0.2 pc, a value ^ 3 times larger than the 
tidal radius estimated by Bergond, Leon & Guibert iiwni 
and Tadross et al. We note that in the present 

work the tidal radius could not be calculated since the 
nonlinear least-squares routine did not reach convergence 
with the three-parameter King model applied to 

the background-subtracted radial distribution of stars in 
NGC 2516 (Fig.El). 

3.2.2. Mass functions 

Similarly to NGC 2477 we derive MFs for the following 
internal regions of NGC 2516: (i) 0.0 < r(pc) < 0.6 (core), 
(ii) 0.6 < r(pc) < 1.7, (Hi) 1.7 < r(pc) < 6.2 and (iv) 
0.0 < r(pc) < 6.2 (overall). The colour-magnitude filter 
(Fig.El) has been applied to the CMDs of the above regions 
and offset held (7.6 < r(pc) < 10.18) to take into account 
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Fig. 4. Left panel: J x (J — H) CMD of the central 20' 
(~ 2.5 pc) of NGC 2516 with the best-fit Padova isochrone 
(age ~ 160 Myr) superimposed. Right panel: CMD of 
the corresponding (same area) offset field. The colour- 
magnitude filter is shown in both panels. Note that the 
MS of NGC 2516 is less populated than that of NGC 2477 
in Fig^ 

the Galactic contamination. The resulting MFs covering 
the mass range 0.39 < m(M 0 ) < 4.17 are shown in Fig.|Hl 

Parameters derived from the MFs, such as observed 
and extrapolated mass, projected and mass densities in 
each region of NGC 2516 are given in Table |21 Similarly 
to NGC 2477 a break in the MFs can be seen, particu¬ 
larly in the core region at m = 0.9 Mq. In addition, mass- 
segregation effects are reflected on the slope variation with 
distance to the center of NGC 2516. In the mass range 
0.90 < m(MQ) < 4.17 the MF slope is similar (y 0.7) 
both for the core and the region 0.6 < r(pc) < 1.7, increas¬ 
ing to X « 1.8 for 1.7 < r(pc) < 6.2. In this mass range 
the overall MF is similar to a standard Salpeter law. In 
the range 0.39 < m(M 0 ) < 0.90 the core MF is very flat 
(x « ~ 1 - 1 ) while in the other regions the slope is similar 
to that of the overall MF (x ~ 0.4). 

The sharp flattening in the MFs of NGC 2516 was pre¬ 
viously detected by Jeffries, Thurston & Hambly mm 
based on CCD photometry of the spatial region internal to 
r « 31' (~ 4pc). In the range 0.3 < m(M 0 ) < 0.7 they de¬ 
rive MF slopes of X = —0.75±0.20or x = —1.00±0.18, de¬ 



0 1 10 


R{pc) 

Fig. 5. Same as Fig. El for the radial density profile of 
NGC 2516. 

pending on the (solar metallicity) isochrone models used. 
In the range 0.7 < m(M 0 ) < 3.0 they derive slopes 
X = 1.47 ± 0.11 or X = 1-58 ± 0.11, for the same models 
as above. Considering the slope uncertainties and differ¬ 
ences in the spatial regions sampled, the MF slopes derived 
in the present work are consistent with those of Jeffries, 
Thurston & Hambly mm . in both mass ranges. In ad¬ 
dition, the total mass derived in the present work m^ot = 
1280 ± 162 M 0 agrees with the value estimated by Jeffries, 
Thurston & Hambly II2001II . m « 1240 — 1560M 0 (cor¬ 
rected for mass segregation for stars with m > 0.3 M 0 ). 
We note that previous works gave exceedingly low values 
for the stellar mass in NGC 2516, e.g. Bergond, Leon & 
Guibert mm estimate mtot ~ 170 M 0 (for stars with 
m > 1.1 M 0 ), while in Tadross et al. mm ,nitot « 
36 M 0 . Extrapolated and observed core mass agree, while 
the extrapolated overall mass is ~ 40% larger than the 
mass locked up in the observed stars. 

Mass segregation in NGC 2516, as implied by the spa¬ 
tial variation of the MF slopes (Fig.E)), is consistent with 
the overall t^eiax 90 Myr and treiax 0.4 Myr in the core, 
both estimates are smaller than the age of the cluster. 

3.3. Comparing dynamical states 

The less-advanced dynamical state of NGC 2516 as com¬ 
pared to that of NGC 2477 is reflected both in the ab- 
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Fig. 6. Mass functions in different spatial regions of 
NGC2516. Each panel contains MFs derived indepen¬ 
dently from the J, H and Ks 2MASS photometry. A break 
in the MFs at mureak ~ 0.9 M© can be seen particularly 
in the core region (panel (d)). ME fits {4>{m) oc 
are shown as dashed lines. 



m(MJ 


solute value of the ME slopes and in the spatial rate of 
change of y. In the mass range 0.90 < m(MQ) < 4.17 of 
NGC2516 the ME slope varies from y 0.7 in the core 
to X « 1.8 in the region 1.7 < r(pc) < 6.2, while in the 
equivalent mass range 1.08 < m(M 0 ) < 1.95 of NGC 2477 
the slope varies from the flatter x ~ 0.0 in the core to 
the steeper x ~ 3.3 in the region at 5.2 < r(pc) < 11.6. 
Thus, in this mass range the radial variation of the ME 
slope in NGG2516 is ^ = (0.30 ± 0.19) pc“^, while in 
NGG2477 = (0.43 ± 0.19) pc“^. A similar situation 

with respect to ME slope values and radial variations oc¬ 
curs in the mass range corresponding to the ME breaks 
in NGC 2516 with ^ = (0.4 ± 0.3) pc"! and NC2477 
with ^ = (0.5 ± 0.2) pc“^. Finally, the overall ME slopes 
in NGC 2516 and NGC 2477 are essentially the same in 
the low-mass MS (x ~ 0.4) and in the high-mass MS 
(X « 1.4). 

Thus, the discussions above show that the fact that 
NGC 2477 is ^ 7 times as old as NGC 2516 turns out to be 
more important to the internal dynamical evolution than 
NGC 2477 being ^ 4 times more massive than NGC 2516. 


4. Multi-parametric analysis 

The analysis method described in Sects. 13.11 and 13.21 re¬ 
spectively for NGC 2477 and NGC 2516 was applied to 
the remaining clusters in Tabled We derived a homoge¬ 
neous set of parameters associated to the internal struc¬ 
ture (core and limiting radii), stellar evolution (extrap¬ 
olated mass) and dynamical evolution (core and overall 
mass and mass density, and MF slope). We do not in¬ 
clude in the statistical analyses below the tidal radius, 
since convergence in the nonlinear least-squares routine 
with the three-parameter King model has been 

reached only for the background-subtracted radial pro¬ 
files of M26 (Rtidai = 15 ± 7pc - 1.9 X Rii„i), NGC 2477 
(Rtidal = 16±2pc^ 1.4xRiim), M67 (Rtidai = 16±3pc^ 
1.4 X Rlim) and NGC 188 (Rtidai = 22 ±3 pc - 1.8 x Ru^). 
As a comparison between conceptually different estimates 
of cluster dimension, in these four open clusters the ra¬ 
tio Rtidai/Riim varies from 1.4 to 1.9. Thus, for sample 
completeness purposes we use the limiting radius in what 
follows. 

Except for cluster age and Galactocentric distance 
which are listed in Tabled the remaining parameters are 
listed in Table 121 in which we consider separately the core 
and overall regions. By columns: (1) - cluster designation; 
(2) - core radius; (3) - observed core mass, which includes 
the mass of the evolved stars down to the mass correspond¬ 
ing to the 2MASS faint magnitude limits ISect. rTL^ : (4) - 
extrapolated core mass, including stars with m = 0.08 Mq; 
(5) - projected core mass density (m/7rR^Q,,g); (6) - core 
mass density (m/|7rRpQj.g); (7) - core MF slope (for the 
mass range m > mbreak); (8) - limiting radius; (9) - 
observed overall mass; (10) - extrapolated overall mass; 
(11) - overall projected mass density; and (12) - over¬ 
all mass density. The overall MF slope (m > mbreak) is 
given in col. 13. In col. 14 we give mbreak) the stellar mass 
where the MF break occurs. Finally, in col. 15 we give 
c* = log (Riiin/Rcore), which Can be measured essentially 
in all open clusters. Derived from the limiting radius, c* 
is not the same as the concentration parameter of glob¬ 
ular clusters, which is based on the tidal radius (Trager, 
King & Diorgovski ll995ll . Since the tidal radius can only 
be measured for the most populous open clusters and/or 
those at high Galactic latitudes (see above), we empha¬ 
size the importance of using in the future c* for statistical 
purposes when dealing with samples of open clusters. The 
relatively narrow range of c* from « 0.8 to ~ 1.0 is char¬ 
acteristic of loose star clusters (Trager, King & Djorgovski 

CMi. 

To characterize the sample we present average values of 
some parameters. For the massive (m > 1 000 Mq) clusters 
we obtain (Rcore) = 1.0±0.3pc, (Mgore) = (3±2 )x10^Mq, 
(Xcore) = 0.5 ±0.5, (Rii„,) = 8.8 ±2.3 pc, (Moverall) = 
(2.8 ± 1.5) X 10^ Mq and (xoveraii) = 1-5 ± 0.2. For the 
less-massive clusters they are (Rcore) = 0.9 ± 0.2 pc, 
(Mcore) = (0.8 ±0.4) X 102 Mq, (Xcore) = -0.9 ±1.1, 
(Rlim) = 7.9±2.5pc, (Moveraii) = (0.7±0.2) X 103 Mq and 
(Xoveraii) = 0.7±0.7. On average the core and overall MFs 
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Table 3. Structural and dynamical-related parameters of the target open clusters. 


Cluster 

(1) 


Core 





Overall 



mbreak 

(Mq) 

(14) 

c* 

(15) 

Rcore ^^obs 

(pc) (lO^Mo) 

(2) (3) 

mtot 

(lO^Mo) 

(4) 

(7 

(5) 

P 

(6) 

X 

(T) 

Riim 

(pc) 

(8) 

mobs 

(IO^^Mq) 

(9) 

mtot 

(ib'Mo) 

(10) 

<7 

(11) 

P 

(12) 

X 

(13) 

M26 

0.8 ± 0.1 1.2 ± 0.3 2.9 ± 0.9 

150 ± 48 142 ± 46 

1.3 ± 0.2 

8.0 ± 0.5 

6.8 ± 1.6 

19 ± 7 

9.4 ± 3.3 

0.9 ± 0.3 

1.5 ± 0.2 

(t) 

1.0 

NGC 2516 

0.6 ± 0.1 0.9 ± 0.1 

0.9 ± 0.2 

82 ± 14 

103 ± 17 

0.7 ± 0.2 

6.2 ± 0.2 

9.1 ± 0.8 

13 ± 2 

10.5 ± 1.3 

1.3 ± 0.2 

1.4 ± 0.1 

0.90 

1.0 

NGC 2287 

1.1 ± 0.1 0.8 ± 0.2 0.9 ± 0.2 

22 ± 5 

15 ±4 

0.2 ± 0.3 

7.0 ± 0.7 

5.7 ± 0.7 

7.2 ± 0.9 

4.7 ± 0.6 

0.5 ± 0.1 

1.5 ± 0.2 

1.01 

0.8 

M48 

0.9 ± 0.2 0.5 ± 0.1 

0.6 ± 0.1 

26 ± 3 

23 ± 3 

0.3 ± 0.2 

8.8 ± 0.4 

4.8 ± 0.7 

6.3 ± 1.5 

2.6 ± 0.6 

0.22 ± 0.05 

1.3 ± 0.2 

1.06 

1.0 

M93 

0.7± 0.1 1.1 ± 0.1 

1.4 ± 0.1 

85 ± 8 

88 ± 9 

-0.2 ± 0.2 

7.2 ± 0.6 

7.4 ± 0.6 

17 ± 5 

10.2 ± 3.4 

1.0 ± 0.3 

1.2 ± 0.1 

— 

1.0 

NGC 5822 

1.1 ± 0.1 1.1 ± 0.1 

1.6 ± 0.4 

43 ± 11 

30 ± 7 

0.6 ± 0.2 

8.0 ± 0.4 

13.6 ± 1.0 

26 ± 10 

13 ± 5 

1.2 ± 0.4 

1.7 ± 0.1 

— 

0.9 

NGC 2477 

1.4 ± 0.1 4.1 ± 0.4 5.2 ± 0.8 

85 ± 13 

45 ± 7 

0.0 ± 0.2 

11.6 ± 0.7 

27 ± 14 

53 ± 16 

13 ± 4 

0.8 ± 0.2 

1.5 ± 0.1 

1.08 

0.9 

NGC 3680 

0.5 ± 0.2 0.3 ± 0.1 

0.4 ±0.1 

45 ± 6 

64 ± 9 

-1.1 ± 0.4 

5.8 ± 0.6 

1.7 ± 0.2 

3.6 ± 1.1 

3.4 ± 1.0 

0.5 ± 0.1 

0.5 ± 0.3 

— 

1.0 

IC4651 

0.8 ± 0.1 0.9 ± 0.1 

1.0 ± 0.1 

54 ± 6 

53 ± 6 

-1.4 ± 0.3 

6.0 ± 0.3 

7.3 ± 0.7 

8.4 ± 1.0 

7.5 ± 0.9 

0.9 ± 0.1 

0.6 ± 0.3 

0.97 

0.9 

M67 

1.1 ± 0.1 1.3 ± 0.2 

1.4 ± 0.2 

34 ± 4 

22 ± 3 

-2.5 ± 0.4 

11.7 ± 0.6 

9.4 ± 0.7 

9.9 ± 1.2 

2.5 ± 0.3 

0.16 ± 0.02 

-0.3 ± 0.2 

0.80 

1.0 

NGC 188 

1.3±0.1 0.7±0.1 

6.1 ± 0.3 

99 ±42 

53 ± 23 

0.6 ± 0.7 

11.9 ± 0.5 

3.8 ± 0.1 

38 ± 16 

8.5 ± 3.6 

0.54 ± 0.23 

1.9 ± 0.7 

(t) 

1.0 


Table Notes, mobs mtot in cols. 3, 4, 9 and 10 includes the mass in the evolved stars mass, mbreak in col. 14 is the mass where the MF break 
begins, x in cols. 7 and 13 is the MF slope in the mass range m > mbreak- (t) ■ because of the distance from the Sun, the low-mass end of the 
MF of M 26 occurs at m = 0.98 Mq. (|) - because of the distance from the Sun and old age the MF of NGC 188 is restricted to the mass range 
0.9 Mq — I.OMq. The uncertainty in mbreak in col. 14 is ~ 0.05 Mq in all clusters. The uncertainty in c* in col. 15 is 0.1 in all clusters except 
for NGC 3680 for which it is 0.2. Units of <7 (cols. 5 and 11) and p (cols. 6 and 12) are Mq pc~^ and Mq pc“^, respectively. 


of the less-massive clusters tend to be flatter than those of 
the massive ones. Core and overall masses of the massive 
clusters are ^ 4 times larger than those of the less-massive 
ones. Within the standard deviations core and limiting 
radii are similar on average, in both types of clusters. 

Because of the reduced number of clusters in this work, 
the results below should be taken as a sampling of the class 
of open clusters. On the other hand, the statistical signif¬ 
icance and homogeneous method from which the parame¬ 
ters were derived guarantee that the results may serve as 
indicative of trends and correlations. 

4.1. Cluster structural and evolutionary parameters vs. 
cluster age 

Relations involving cluster parameters and cluster age are 
investigated in Fig. 0 in which we use different symbols 
to identify massive and less-massive clusters. Among the 
parameters tested the best correlation with age occurs for 
Riim (panel (a)). This correlation should be expected since 
young (age < 700 Myr, Friel 1199^1) clusters with large di¬ 
ameters probably are not bound systems (Lynga 119821 
Janes, Tilley & Lynga EHHl), and small, massive clusters 
will be dissolved due to the strong effects of the internal 
dynamics (Tadross et al. 1200211 . This correlation was not 
present in the sample studied by Tadross et al. (Eoni. 
The core radius in panel (g) does not seem to correlate 
with age. 

An interesting pattern occurs in the plots Pcore x age 
(panel (1)) and acore x age (panel (i)) which may reflect dif¬ 
ferent evolutionary processes in the core of more-massive 
clusters with respect to less-massive ones. A similar pat¬ 
tern occurs in the plot Xoveraii X age (panel (d)), in which 
massive and less-massive clusters seem to evolve differ¬ 
ently, in the sense that old, less-massive clusters tend to 
have flatter overall MFs than old massive ones. Because 


of their looser nature with respect to massive clusters, 
old less-massive clusters have lost a significant fraction 
of the low-mass content through internal processes such 
as mass segregation and evaporation, and external ones 
such as tidal stripping, disk shocking and encounters with 
molecular clouds (Bergond, Leon & Guibert 1200111 . As a 
consequence, a higher flattening degree in the MFs of old 
less-massive clusters with respect to old massive ones, is 
expected. A similar pattern occurs with respect to the re¬ 
lation of the core MF slope with age in panel (j). 

The projected overall mass density (panel (c)), total 
core mass (panel (k)) and overall mass density panel (f)) 
do not correlate with age, but massive and less-massive 
clusters in general occupy distinct positions in the plots. 

Considering the significant scatter in the plots we find 
no indication of correlation with cluster age of total overall 
mass (panel (e)), observed overall mass (panel (b)) and 
observed core mass (panel (h)). A larger cluster sample is 
necessary to check these points. 

4.2. Cluster structural and evolutionary parameters vs. 

Galactocentric distance 

Relations involving cluster parameters and Galactocentric 
distance are investigated in Fig. 0 The plot in panel (f) 
suggests that there is no correlation of cluster mass with 
Galactocentric distance. Although the scatter in most 
plots is considerable, the plot in panel (e) suggests a cor¬ 
relation of Riim with Galactocentric distance dec j similar 
to the results of Lynga and Tadross et al. mm . 

This correlation holds for massive and less-massive clus¬ 
ters. The core mass density (panel (a)) and, to a less extent 
the overall mass density (panel (i)), seem to fall for clus¬ 
ters at larger deC) probably because of the slight depen¬ 
dence of core radius on dec (panel (g)). Finally, in panel 
(i) of Fig. 0 we examine the relation of dec with cluster 
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Fig. 7. Relations of cluster parameters with cluster age. 
Filled symbols correspond to massive clusters (overall 
mass > 1 000 Mq). 


age. Young and old clusters, with any mass, are similarly 
distributed in terms of doc- This apparently contradicts 
the conclusion of Lynga in which older clusters are 

found preferentially in the outer parts of the Galaxy, while 
younger clusters are evenly distributed. However, our re¬ 
sult is based on 11 clusters only, which sample a limited 
range in dec in the solar surroundings. 

4.3. Core vs. overall parameters 

In Fig. 13 we examine relations involving core and over¬ 
all parameters. The tight correlations involving core and 
limiting radii (panel (g)) and core and overall total mass 
(panel (d)) are probably consequences of size and mass 
scaling, in the sense that bigger clusters tend to have big¬ 
ger cores and larger masses. Linear least-squares fits to the 
points resulted in correlation coefficients of 0.85 and 0.84, 
respectively for the plots in panels (d) and (g). The pro¬ 
jected overall density (panel (a)) and overall mass density 
(panel (h)) seem to correlate as well, although with sig¬ 
nificant scatter. Core radius and overall mass are tightly 
correlated (panel (f)) as well. However, massive and less- 
massive clusters seem to follow separate, but nearly par¬ 
allel, paths in the plot. Correlation coefficients are 0.97 


3 



dGc(kpc) dQc(kpc) dQc{kpc) 

Fig. 8. Relations of cluster parameters with 
Galactocentric distance. Symbols as in Fig. [3 

and 0.68, respectively for the hts of the massive and less- 
massive clusters. 

Total overall mass and overall MF slope correlate 
(panel (c)). In this case the MS of massive clusters is char¬ 
acterized by a narrow range of MF slopes, 1 < y < 2, while 
the MS of less-massive clusters has y spread over a large 
range. No definite pattern emerges from the relations of 
total overall mass with core mass density (panel (b)) and 
core projected density (panel (e)). 

Finally, the core and overall MF slopes are tightly cor¬ 
related (panel (j)), with a correlation coefficient of 0.90. 
Total core mass and core MF slope (panel (k)) present a 
similar pattern as total overall mass and overall MF slope 
(panel (c)). 

5. MF break and the universality of the IMF 

Since the early work of Salpeter (IM3) in which the 
IMF of solar-neighbourhood stars in the range 0.4Mq to 
IOMq was described by the function tpim) oc 
with y = 1.35, evidence has been accumulating pointing 
to the presence of breaks in the MF of low-mass stars, 
in the sense that it tends to be flatter than the MF of 
high-mass stars. On the high-mass end, the Salpeter IMF 
has been shown to extend up to very-massive stars, such 
as m 120 Mq (Scalo 119861 Massey IM- However, 
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Fig. 9. Relations involving core and overall parameters. 
Symbols as in Fig. [3 Dashed lines correspond to linear 
least-squares fits to the points. 


star counts in the solar neighbourhood show that the 
IMF for the low-mass range breaks up in two power-laws, 
X= 1.2±0.3for 0.5 < m(MQ) < 1.0 and y = 0.3± 0.7 for 
0.08 < m(M0) < 0.5. In the brown dwarf mass range the 
slope seems to be even flatter, y « — 0.5 CKroupa 120011 

Present-day data suggest that the IMF may have a 
universal character, in the sense that it probably is sim¬ 
ilar in very different environments, from star-forming re¬ 
gions in molecular clouds to rich star clusters and the field 
(Kroupa 1200211 . As additional observational evidence of 
similar IMFs in different environments von Hippel et al. 
cnna discuss star counts in the near-IR, in the optical 
with Hubble Space Telescope data and from the Hubble 
Deep Field which indicate that field stars of the Galactic 
thick disk and halo have a LF similar to that of the solar 
neighbourhood down to 0.08 Mq stars - stellar populations 
differing significantly in local density and chemical abun¬ 
dances. As pointed out by Kroupa Ilyina this scenario 
appears to indicate that the distribution of stellar masses 
should depend only on the process of molecular cloud frag¬ 
mentation. In this case, the fragmentation would have to 
produce similar IMFs despite very different initial condi¬ 
tions, a physical process which is not yet fully understood 
(Kroupa r2()()2ll . 


Fig. 10. Open clusters with a MF break in the low-mass 
MS. Because of the significant error bars in the MF of 
NGC 2287 (panel (b)) the low-mass MS and high-mass MS 
fits do not match at m « 1 Mq. Gore and overall regions 
are considered separately. 

The relatively populous nature of the open clusters 
included in the present study offers an opportunity to in¬ 
vestigate the presence and properties of the MF break 
in clusters with a high statistical significance. However, 
the basically solar metallicity of the present clusters does 
not allow inferences on relations of mbreak with metallicity 
(e.g. von Hippel et al. 11991)11 or on the possible dependence 
of MF flattening (in the range m < 0.7 Mq) with metal¬ 
licity, as discussed in Kroupa II2()()2II . 

An MF slope break, followed by a MF slope flatten¬ 
ing (turnover) for lower masses, is present in six clus¬ 
ters of the present sample. The MFs of NGG 2477 clearly 
present a break at mbreak ~ 1.08 Mq (Fig. |3), while in 
NGG 2516 the break occurs at mbreak ~ 0.90 Mq (Fig. EJ. 
In both clusters the break consistently occurs in the four 
spatial MFs presented in Figs. 01 and 01 In the remain¬ 
ing clusters of the present paper the break is present 
in the MFs of M48 at mbreak ~ 1.06 Mq (overall MF), 
NGG 2287 at mbreak ~ 1.01 Mq (core and overall MF), 
IC4651 at mbreak ~ 0.97 Mq (overall MF) and M67 at 
mbreak ~ 0.80 Mq (core and overall MF). The core and 
overall MFs of these clusters are shown in Fig. HOI along 
with the respective MF fits. In col. 2 of Table 0) we give 
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Table 4. Core and overall MF slopes for m < nibreak- 


Cluster 

Am 

X(core) 

X (overall) 

(1) 

(Mo) 

(2) 

(3) 

(4) 

NGC 2516 

0.40 - 0.90 

-1.1 ±0.5 

±0.4 ±0.2 

NGC 2287 

0.56 - 1.01 

-1.3 ±0.5 

-0.2 ±0.3 

M48 

0.60 - 1.06 

±0.3 ±0.2 

0.0 ±0.7 

NGC 2477 

0.76 - 1.08 

-2.3 ±0.5 

±0.4 ±0.4 

IC4651 

0.63 - 0.97 

-1.4 ±0.3 

-1.4 ±0.7 

M67 

0.66 - 0.80 

-3.9 ± 1.1 

-2.5 ±0.2 

Table Notes. 

Col. 2: Mass 

range where the 

MF flatten- 


ing is observed. Cols. 3 and 4: Mass-function slope 
(<?!)(m) oc for m < mureak- 


the mass range over which we measure the MF flattening 
(m < mbreak), while the corresponding core and overall 
MF slopes in this mass range are in cols. 3 and 4, respec¬ 
tively. 

The flattening degree in the MF for masses below 
mbreak IS indeed quite high (Table 0)) when compared to 
that in the mass range above mbreak (Table 13, particu¬ 
larly in the core. Mass segregation throughout the cluster 
and low-mass star loss through dynamical evaporation cer¬ 
tainly contribute to MF flattening. This effect is expected 
to be more significant in the old and less-massive clus¬ 
ters (Bergond, Leon & Guibert [200111 . In the mass range 
m > mbreak this effect is confirmed by the relations y x age 
both in the core and overall MFs (Fig. [7| panels (j) and 
(d), respectively). In addition, observational consequences 
of the presence of a significant fraction of undetected bina¬ 
ries (and higher-order multiple systems) in a cluster are, 
among others, the underestimation of the number of low- 
mass with respect to higher-mass stars and widening of 
the MS. Consequently, somewhat flat MFs are expected 
to occur in clusters with a high fraction of unresolved 
binaries (e.g. Bonatto, Bica & Santos Jr. and ref¬ 
erences therein). Indeed, based on proper-motion consid¬ 
erations Bica & Bonatto derived that the fraction 

of undetected binaries in the cores of M26, NGC2287, 
M48, M93, NGC5822, NGC3680, IC4651 and M67 is 
in the range 15% to 50%. In the core of the old open 
cluster NGC 188 Bonatto, Bica & Santos Jr. (120(1411 es¬ 
timated that the fraction of undetected binaries is prob¬ 
ably around 100%. However, even fractions of ^ 100% 
of unresolved binaries cannot produce such flat MFs as 
those detected above (Bonatto, Bica & Santos Jr. 12004(1 . 
Thus, the flat MFs and the presence of the MF break in 
clusters of any age and mass (Table 13 and Fig. must 
essentially reflect the effects of the internal dynamics of 
clusters on the MFs and/or some fundamental property 
of the IMF associated to different conditions in star for¬ 
mation. Kroupa points out that systematic differ¬ 

ences in the IMF should appear between low-density en¬ 
vironments and high-density regions. To investigate this 
issue we test relations involving parameters of the clusters 
which present the MF break. The results are in Fig. 
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Fig. 11. Relations involving parameters of the clusters 
with the MF break. Xcore and Xoveraii in panels (a), (d), 
(e), (g) and (h) are the MF slopes derived for the range 
m < nibreak- Symbols asin Fig. [T] 


Systematic MF flattening with time due to dynamical 
evolution is manifest in panel (a) of Fig. 1111 which shows 
that the overall MF slope (in the mass range m < mbreak) 
decreases with age, particularly for the less-massive clus¬ 
ters. A similar trend is also present in the core MF slope 
(panel (d)). The decrease of core and overall MF slopes 
with age was present in the mass range m > mbreak (Fig.EI 
panels (j) and (d)). 

Similarly to the mass range m > mbreak) core and over¬ 
all MF slopes correlate as well (panel (g)), particularly in 
the less-massive clusters. The scatter hinders conclusions 
on dependences of MF slope on mass density in the core 
(panel (h)) and in the whole cluster (panel (e)). Except for 
the populous cluster NGC 2477, mbreak seems to decrease 
with overall cluster mass (panel (b)). The same trend oc¬ 
curs for core mass, since core and overall mass are tightly 
correlated (panel (j) in Fig. 0. Within the scatter we do 
not find correlations of mbreak with mass density, both in 
the core (panel (f)) and in the whole cluster (panel (c)). 
However, mbreak seems to decrease with cluster age for the 
clusters older than 1 Gyr (panel (i)). 

MF slopes in the mass range m < mbreak (break region) 
and m > mbreak correlate, both in the core (panel (j)) and 
particularly in the whole cluster (panel (k), correlation 
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Fig. 12. Open clusters with no clear evidence of a MF 
break. Because of the distance from the Sun, the MF of 
M26 in panel (a) reaches down to m « 0.98M©. 

coefficient 0.96). This suggests a common origin for the 
mechanisms responsible for the excess MF flattening in 
m < mbreak and that in m > mureak- 

In Fig. El we show the MFs of the clusters which 
can be fitted by a single power-law. Because of the rel¬ 
atively large distance from the Sun (Table 0 , the low- 
mass end of the observed MF of M26 (panel (a)) oc¬ 
curs at m « 0.98 Mq, nearly the same value in which 
the MF break usually begins. Deeper observations are 
needed to access the MF break region in M 26. For M 93 
(panel (b)) and NGC3680 (panel (d)) the single power- 
law extends down to « 0.7 Mq and ~ 0.8 Mq respec¬ 
tively, with no evidence of a break. This is particularly 
true for the MF of NGC5822 (panel (c)), which extends 
to « 0.53 Mq fitted with a single power-law. The MF of 
NGG 188 is not shown in Fig. El because the distance 
from the Sun (dQ ~ 1.7 kpc) combined to the old age 
(« 7 Gyr) produce a MF restricted to the narrow mass 
range 0.9 < m(MQ) < 1.1. The mass range is too narrow 
for any inferences on the presence of a MF break in this 
old open cluster. 

Recall that the MFs of the present cluster sample 
have been built by taking into account the 2MASS 99.9% 
Completeness Limit ISect. Consequently, the low- 

mass end of each MF presented in this paper corresponds 


to the same apparent magnitudes, J = 15.8, H = 15.1 and 
Ks = 14.3. In this sense, if the MF break was artificially 
created by detection incompleteness in the faint magni¬ 
tude limit, it should be present in the MFs of all clusters, 
irrespective of the corresponding low-mass limit. This is 
not the case of the MFs of M26, M93, NGG 5822 and 
NGC3680 (Fig. El), in which the low-mass limit varies 
from « 0.5 Mq (NGC 5822) to « 1.0 Mq (M26). 

he presence of the MF break is not associated to clus¬ 
ter age and mass or concentration parameter, at least in 
the studied clusters. However, the value of mbreak seems 
to decrease with cluster mass. We conclude that the MF 
break for m < 1 Mq is definitively present in NGC 2516, 
NGC 2287, M48, NGC 2477, IC4651 and M67. On the 
other hand it is absent, at least for the mass range 
m > 0.7 Mq, in M93, NGC 5822 and NGC 3680. Further 
research accessing low-mass stars in a statistically signif¬ 
icant sample of open clusters is necessary to settle the 
issue whether or not the MF break is present in all cases, 
and thus test the hypothesis of the universal character of 
Kroupa’s IMF. 

5.1. Mass segregation and low-mass stars evaporation 

A comparison of the core and overall MF slopes in Table 131 
may help us to evaluate the amount of mass segregation 
and low-mass stars evaporation which have already took 
place in the clusters. Except for M 26, the core MF slope 
in the remaining clusters is much flatter than the overall 
one. In the young and rather massive (Table |31) cluster 
M26 both MF slopes correspond essentially to the stan¬ 
dard Salpeter one, within uncertainties. This suggests that 
mass segregation and low-mass stars evaporation have not 
yet had time to produce observable changes neither in the 
internal structure of M 26 nor in the MFs. This fact is con¬ 
sistent with the overall relaxation time treiax 180 Myr, 
which is larger than the cluster age (^ 70Myr). On the 
other hand, the similarity between the core and overall 
MF slopes in M26 occurs despite the much shorter core 
relaxation time of treiax 3.2 Myr. This apparent contra¬ 
diction reflects the complex way in which internal dynam¬ 
ical evolution relates to age and relaxation time, and in 
particular the different evolution time-scales associated to 
the core and overall regions. 

In order to investigate the relationship between relax¬ 
ation time and cluster age with dynamical evolution, and 
to estimate the corresponding effects on MFs we calculate 
for each cluster the evolutionary parameter (t), which is 
defined as the ratio of the cluster age to the relaxation 
time, T = age/treiax- We take advantage of the 2MASS 
spatial coverage to examine separately the evolution of 
the core and overall regions. In Table [31 we list treiax and 
r for each cluster. 

Considering separately massive and less-massive clus¬ 
ters, we found that core and overall relaxation times are 
linearly related. For the massive clusters we derive the re¬ 
lation treiax (overall) ~ (89 ± 15) x treiax(core) (correlation 
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Table 5. Relaxation time and evolutionary parameter. 


Cluster 

(1) 

Core 

Overall 

treiax 

(Myr) 

(2) 

T 

(3) 

treiax 

(Myr) 

(4) 

T 

(5) 

M26 

3.2 ±2.0 

22± 14 

178 ± 116 

0.4 ±0.3 

NGC 2516 

0.4 ±0.1 

365 ± 74 

94± 24 

1.7 ±0.4 

NGC 2287 

0.8 ±0.2 

212 ±46 

40± 10 

4.0 ± 1.0 

M48 

0.7 ±0.2 

517 ± 166 

48 ±28 

7.0 ±4.5 

M93 

0.9 ±0.2 

457±112 

147 ± 93 

2.7 ± 1.7 

NGC 5822 

2.4 ± 1.3 

423 ± 234 

262 ± 170 

3.8 ±2.5 

NGC 2477 

6.2 ±2.5 

181 ± 75 

597 ± 344 

1.9 ± 1.0 

NGC 3680 

0.2 ±0.1 

7100 ± 2700 

29± 17 

55 ±32 

IC 4651 

0.6 ±0.1 

2854 ± 455 

33 ±7 

54± 13 

M67 

1.2 ±0.2 

2611 ± 402 

80± 15 

40 ±8 

NGC 188 

14 ±9 

504 ± 321 

631 ± 400 

11 ± 7 


Table Notes. The relaxation time is calculated according to 
t.eiax « 0.12 (j^) (^) (Myr), where N is 


the total number of stars, R is the radius of the region 
considered and is the velocity dispersion. We adopted 
(Tv ~ 3kms“^ (Binney & Merrifield I1998t . r = age/t^eiax 
is the evolutionary parameter. 


coefficient of CC ~ 0.94) and an average of ~ 120. For the 
less-massive ones treiax(overall) ^ (56 ± 20) x treiax(core) 
(CC « 0.85) and an average of ^ 80. 

We assume that flat overall MF slopes can be ac¬ 
counted for by low-mass stars evaporation resulting from 
large-scale mass segregation and external effects such as 
tidal stripping by the Galactic gravitational field. In this 
case, overall MF flattening is expected to be associated 
to r(overall). In fact, in the clusters M26, NGC2516, 
NGC2287, M48, M93, NGC5822 and NGC2477, which 
have overall MF slopes similar to that of Salpeter 
(thus, little low-mass stars evaporation), r(overall) < 7. 
Consequently, the age of these clusters results smaller than 
~ 7 X treiax (overall) which turns out to be too short a time 
for the evaporation to produce significant changes in the 
overall MFs. In this sense the flat core MFs of NGC 2516, 
NGC2287, M48, M93, NGC 5822 and NGC 2477 can be 
accounted for by mass segregation, since for these clusters 
r(core) > 200. The nearly Salpeter MF in the core of M 26 
is probably due to the small r(core) ^ 22. 

On the other hand, NGC 3680, 1C 4651 and M67, 
which have flat overall MF slopes, have T(overall) > 40. 
The very flat core MFs in these clusters result from 
r(core) > 2600. Thus, the core and overall evolutionary 
parameters are large enough for the large-scale mass seg¬ 
regation and low-mass stars evaporation to have produced 
significant MF flattening in these clusters. 

Although with r (overall) ^11 and T(core) ^ 504, the 
oldest cluster in this work NGC 188 doesn’t seem to fit in 
the above scenario, because its overall MF slope is x 1-9. 
However, because of the old age and large distance from 
the Sun we could access the MF only over the reduced 




Fig. 13. MF slopes as a function of the evolutionary pa¬ 
rameter T. Left panels: MFs in the range m > mbreak; right 
panels: MFs in the range m < mbreak- Top panels: the scale 
in the abscissa is based on the overall r; bottom panels: 
abscissa in terms of core r. Dashed lines in panels (a), (b) 
and (c): tentative fit with x('r) = Xo ~ Xi exp (— tq/t). 
Symbols as in Fig.[71 

range 0.9 Mq — 1.0 Mq, which produced large uncertainties 
in the MF slope. 

We summarise the evolutionary scenario of core and 
overall MFs with respect to r in Fig. ^1 The overall MF 
slopes (for m > mbreak) of the massive clusters (panel 
(a)) do not change significantly for 0.4 < r(overall) < 11. 
On the other hand, the overall MFs of the less-massive 
clusters begin to flatten for r(overall) > 7, which means 
that evaporation effects tend to become appreciable on 
the MFs when the age of the cluster becomes larger than 
^ 7 X treiax (overall). This occurs probably because of the 
looser nature of the less-massive clusters. 

Because of high density and small dimensions (and 
thus small treiax), core MFs tend to flatten because of mass 
segregation for r(core) > 22 (M26) and T(core) < 150, in 
massive and less-massive clusters alike (panel (c)). To first 
order we take T(core) = 100 as representative of the value 
for which core MF flattening becomes significant. 

Based on the discussion above we can express the time 
it takes for significant changes in the MFs (in the range 
ni > mbreak) to be produced as At = r x treiax- In the less- 
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massive clusters the core and overall MFs begin to flatten 
at r(core) ^ 100 and r(overall) ^ 7, respectively. Taking 
into account the average ratio treiax (overall)/tj-eiax (core) « 
80, we find At(overall) ~ 6 x At(core). We conclude that 
appreciable slope flattenings in the overall MFs of the less- 
massive clusters take on average ~ 6 times longer to occur 
than in the core MFs. 

The value of T(overall) in which the overall MFs 
of massive clusters begin to flatten is not clear in 
panel (a) of Fig. El Thus, assuming it as r(overall) ~ 
11 (NGC 188) for simplicity, and considering the average 
ratio treiax(overall)/treiax(core) « 120, we find that slope 
flattenings in the overall MFs of massive clusters take on 
average ~ 13 times longer to occur than in the core MFs. 

In the mass range m < mureak the overall MFs (panel 
(b)) flatten for r (overall) > 2, which indicates that in 
this mass range dynamical evolution effects are effective 
in flattening the MFs quite rapidly. The scatter in panel 
(d) does not allow to draw conclusions on the behaviour 
of the core MFs. 

The discussion above (and the apparent dependence of 
y on r in Fig. uni suggests that the MF slopes undergo an 
exponential decay with r. Despite the reduced number of 
points we tentatively applied a fit with the empirical func¬ 
tion x(t) = Xo ~ Xi 6xp (—tq/t) to the core and overall y 
distributions (panels (a), (b) and (c) in Fig. EJ, combin¬ 
ing massive and less-massive clusters. For y(overall) we 
derived yo = 1.46 ± 0.10, xi = 1.9 ± 0.7 and tq = 18 ± 14, 
with a correlation coefficient CC = 0.88. For y(core) we 
found yo = 1.17±0.23,yi = 3.0±0.7andTo = 439±156, 
with CC = 0.88. The resulting fits are shown in Fig. El 
For the overall y in the range m < moreak (panel (b) we 
derive yo = 0.42 ± 0.27, yi = 3.4 ± 0.4 and tq = 8 ± 4, 
with CC = 0.98. These analytical relations can be used 
to compare core and overall MF flattening time scales. 
According to the above relations the overall MFs of the 
less-massive clusters would flatten by Ay = 0.35 in ^ 4x 
the time it would take for the core MF to flatten by the 
same amount. In the massive clusters this ratio would be 
~ 6. Both values basically agree with our previous esti¬ 
mates (see above). 

In the same way, the overall MFs in the range m < 
mbreak flatten by Ay = 0.35 in ^ 1/3 of the time the 
m > mbreak Overall MFs flatten by the same amount. This 
is consistent with the fact that on the way towards energy 
equipartition, stars in the range m < mbreak should have, 
on average, a larger velocity dispersion (ctv "^1/-s/m) than 
those in the range m > mbreak and thus, they may leave 
out the cluster first. 

A larger number of clusters is necessary to first check 
the existence of this analytical representation of y(T) and 
then to derive more accurate values of yo, yi and tq, and 
finally, discuss further implications. 

6. Evidence of an open cluster Fundamental Plane 

In Sect. 01 we discussed relations involving pairs of intrin¬ 
sic open cluster parameters. We found a few correlations 


among core and overall parameters, while in some plots 
massive and less-massive clusters are clearly segregated. 
This raises the question whether one of these pairs might 
be related to a third parameter, which could reflect the 
existence of a fundamental plane (FP) of open cluster pa¬ 
rameters, similar to that of the elliptical galaxies (Dressier 
et al. 119871 Djorgovski & Davis EHUl- Elliptical galaxies 
populate a plane in the 3D-parameter space defined by 
the luminosity (L), effective surface brightness (/Je) and 
velocity dispersion (ctv)- The plane comes about because 
as a consequence of the virial theorem a well-defined rela¬ 
tionship between mass (M), /j,e and CTv is expected to occur 
in elliptical galaxies. Scatter in the plane depends essen¬ 
tially on the function relating the ratio M/L with L and /Ue 
(Lucey, Bower & Ellis ESI- There is suggestive evidence 
that the planes of the ellipticals in the Coma cluster and 
that in the Virgo are not parallel (Lucey, Bower & Ellis 



Although complete virialization in open clusters is an 
idealization, as clusters of different masses age they may 
reach advanced dynamical states. Indeed, kinetic theory 
and numerical simulations predict that N-body encounters 
in multi-mass models tend to produce energy equipartition 
within a time-scale of the order of the relaxation time (de 
la Fuente Marcos & de la Fuente Marcos 1200211 . Except 
in M26, the overall relaxation times for the remaining 
clusters in the present work are significantly smaller than 
the cluster age. Core relaxation times are indeed much 
smaller than the ages (Table El . In this sense, varying 
levels of relationship between open cluster parameters are 
expected to occur, particularly in the old clusters. 

Open clusters do differ of elliptical galaxies particu¬ 
larly in terms of stellar population (absolute numbers and 
age), linear dimension and dynamical evolutionary state. 
Thus, the existence of an equivalent FP in open clusters 
raises some questions, ftj Is the open cluster FP related to 
dynamical evolution? fiij Does the massive/less-massive 
cluster segregation (Sect. 01 still hold in the FP? finj 
Does Calactocentric distance affect the FP? A more com¬ 
prehensive sample of open clusters is essential to answer 
these questions. 

In the case of the present open clusters (Tables 01 and 
01 we found that the overall mass and core radius are 
effectively related to the projected overall mass density. In 
Fig.Elwe show a 3D-projection of those three parameters, 
where one can see that they distribute nearly in a plane. 
A similar relation occurs for overall mass, core radius and 
projected core mass density. 

At the present stage the number of open cluster pa¬ 
rameters is of the same order as the number of objects 
already analyzed by the techniques presented in the pre¬ 
vious sections. We are using the same techniques to derive 
parameters of a larger number of clusters with significantly 
different properties, such as age, mass, core and limiting 
radii, densities, Calactocentric distance, etc. Our goal is 
to apply a principal component analysis on this compre¬ 
hensive cluster sample not only to search for correlations 
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Fig. 14. Evidence of a fundamental plane of open cluster 
parameters. The distribution of overall mass (x-axis), core 
radius (y-axis) and projected overall mass density (z-axis) 
essentially defines a plane in this 3D projection. 

among parameters, but to better explore the possibility of 
an open cluster FP as well. 

7. Concluding remarks 

We analysed 11 nearby open clusters with ages in the 
range 70 Myr to 7 Gyr and masses in the range « 400 Mq 
to « 5 3OOM0 based on J, H and Ks 2MASS photom¬ 
etry. The clusters are M26, NGC2516, NGC 2287, M48, 
M 93, NGC 5822, NGC 2477, NGC 3680, IC 4651, M 67 and 
NGC 188. Radial density profiles and mass functions were 
built after taking into account the field contamination. As 
a consequence we derived a homogeneous set of parame¬ 
ters associated to the structure and stellar and dynamical 
evolution of the clusters. 

The method is illustrated in detail by analysing the 
populous open clusters NGC 2477 and NGC 2516 for the 
first time in the near-infrared. For NGC 2477 we derive an 
age of 1.1 ± 0.1 Gyr, distance from the Sun d© = 1.2 ± 
0.1 kpc, core radius Rcore = 1.4 ±0.1 pc, limiting radius 
R-iim = 11.6±0.7pc and total mass (extrapolating the MF 
down to the H-burning mass limit, 0.08 Mq) mtot ~ (5.3± 
1.6) X 10 ^Mq. The MF slope varies significantly in this 
cluster, being flat (% « 0.0) in the core (0.0 < r(pc) < 1.4) 
and steep (y ~ 3.3) in the outskirts (5.2 < r(pc) < 11.7). 
The overall MF has a slope x ~ 1-5. The number-density 
of evolved stars (with respect to MS stars with mass down 
to 0.08 Mq) in the core is about 4 times as large as in 
the cluster as a whole. In the halo the number-density 
of low-mass stars is ~ 4 times as large as that of the 
evolved stars. These facts reflect the effects of large-scale 
mass segregation in NGC 2477. For NGC 2516 we derive 
an age of 160 ± 10 Myr, d© = 0.44 ± 0.02 kpc, Rcore = 
0.6 ± 0.1 pc, Riini = 6.2 ± 0.2 pc and mtot ~ (1-3 ± 0.2) x 
lO^M©. Similarly to NGC 2477, mass-segregation effects 
in NGC 2516 are reflected in the spatial variation of the 
MF slopes, X ~ 0.7 in the core (0.0 < r(pc) < 0.6) and 


X ~ 1.8 in the outskirts (1.6 < r(pc) < 6.2). The overall 
MF has a slope x ~ 1-4. The overall MFs of NGC 2477 
and NGC 2516 have slopes similar to a standard Salpeter 
IMF (x« 1.35). 

Six of the 11 clusters present a break around 1M© fol¬ 
lowed by a sharp flattening in the MF in the mass range 
« 0.5 M© to « IM©. The mass range where we detected 
the MF break in open clusters basically coincides with 
previous results based on star counts in the solar neigh¬ 
bourhood fKroupa 12001112()()2|) . In the clusters analysed 
here the MF break is not associated to mass or concen¬ 
tration parameter. The presence of the MF break in all 
clusters would be consistent with the universal charac¬ 
ter of the IMF ('Krniir)a l2002ll . However, 3 clusters in our 
sample do not present the MF break, at least for the mass 
range m > 0.7 M©. This does not exclude the possibility 
that in some clusters the MF break occurs at lower masses 
or that dynamical effects may somehow damp out features 
in the MFs. To settle these issues we are analyzing a more 
statistically significant sample of open clusters. 

We investigated the dependence of MF slopes on the 
evolutionary parameter r = age/treiax- We found that dy¬ 
namical effects begin to produce substantial flattening in 
the core MFs of massive and less-massive clusters when 
the cluster age is ^ 100 x treiax(core). In the overall MFs 
the changes become noticeable when ^ 7 x treiax (overall) 
and ~ 11 X treiax(overall), respectively for the less-massive 
and massive clusters. Considering the linear relationship 
between core and overall relaxation times, we conclude 
that appreciable slope flattening in the overall MFs of 
the less-massive clusters take ~ 6 times longer to oc¬ 
cur than in the core MFs. In the massive clusters they 
take a time ^ 13 times longer. We found that MF slopes 
vary with t according to the empirical relation xi'’') — 
Xo -Xiexp(-To/r). 

We also searched for relations of cluster parameters 
with cluster age and Galactocentric distance. The main 
results are: (i) cluster size correlates both with cluster age 
and Galactocentric distance; (ii) because of size and mass 
scaling, core and limiting radii, and core and overall mass 
correlate as well; (Hi) core radius and overall mass are 
correlated, but massive (m > lOOOM©) and less-massive 
clusters follow separate, nearly parallel paths on the plane; 
(iv) the MF slopes of massive clusters are restricted to 
a narrow range, while those of the less-massive clusters 
distribute over a wider range; (v) core and overall MF 
slopes are correlated. 

We found relations involving three parameters simul¬ 
taneously which may suggest a fundamental plane of open 
cluster parameters. The relations involve overall mass, 
core radius, and projected overall mass density or pro¬ 
jected core mass density, which are those more closely 
related to the parameters involved in the Fundamental 
Plane of ellipticals. 

As prospective work we are carrying on the present 
analysis for a larger sample in order to draw more quan¬ 
titative conclusions on the fraction of open clusters which 
present a break in the MF, the reason why some do not 
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have the break, and on correlations of mbreak with other 
cluster parameters, further exploring as well the possibil¬ 
ity of a fundamental plane. 
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